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We propose a method for generating and controlling a spatially separated vortex-antivortex pair 
in a Bose-Einstein condensate trapped in a toroidal potential. Our simulations of the time dependent 
Gross-Pitaevskii equation show that in toroidal condensates vortex dynamics are different from the 
dynamics in the homogeneous case. Our numerical results agree well with analytical calculations 
using the image method. Our proposal offers an effective example of coherent generation and control 
of vortex dynamics in atomic condensates. 
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I. INTRODUCTION 

Quantum liquids and gases are an excellent environ- 
ment for studies of solitons and topological quantum 
structures such as vortices. Studies of superfluidity and 
vorticity have been done mostly in Hell (!]]. Due to the 
high density and strong interactions the theoretical stud- 
ies give usually only a qualitative agreement with experi- 
ments. Moreover, the size of the vortex core is very small 
in Hell, typically ~ 1 A. The situation is drastically dif- 
ferent for atomic Bose-Einstein condensates [||. These 
many-body systems are dilute and weakly interacting, 
and as a result a quantitative agreement between theory 
and experiments can often be found. Also, the expected 
size of a vortex core is fairly large, ~ 0.1 — 1 /im, making 
detailed experimental studies possible. 

Several theoretical approaches have been suggested 
for the generation of vortices in condensates, see 
e.g. Refs. PHdQl- Among them, the stirring of a conden- 
sate with a blue-detuned laser in a single component [ pd] | , 
and the coherent interconversion between the two com- 
ponents of a binary condensate [Q have been success- 
fully implemented. In addition, the creation of vortex 
rings after the decay of a dark soliton has been recently 
reported 13 1. These vortices carry one quantum of cir- 
culation each, and they form stable, geometric patterns. 
However, the success in controlling atomic condensates 
suggests that one could create also non-equilibrium situ- 
ations to investigate dynamics of vortices, and in particu- 
lar, vortex collisions. Alternatively, a vortex can be gen- 
erated through the "phase- imprinting" method |Q. This 
approach has been successfully used to generate solitons 
in condensates Q]. 

Here we propose how to create and control a pair of 
spatially separated vortices of opposite circulation via 
controlled decay of dynamically unstable solitons in a 
toroidal trapped condensate. Due to the toroidal geom- 
etry and the employed technique for the creation of the 
vortices, the dynamics and circulation of the generated 
vortices are correlated, even when they are generated at 



well-separated locations. 

In homogeneous condensates and in absence of dissipa- 
tion a vortex-antivortex pair moves as a whole through 
the fluid |Q. Such behavior can be, however, substan- 
tially modified in the case of trapped condensates. In par- 
ticular, for the case of stiff potentials the vortex dynamics 
is altered by the constraint of vanishing normal veloc- 
ity field at the boundaries We show how vortex- 
antivortex dynamics in toroidal geometries become dif- 
ferent from the equivalent dynamics in a homogeneous 
(not trapped) condensate. 

We propose a combination of condensate splitting, 
phase imprinting, and subsequent merging of separated 
parts into a toroidal condensate, which leads to a very 
controlled mechanism of vortex creation. We describe 
the details of our proposal and numerical solution of the 
relevant Gross-Pitaevskii equation in Sec. |h[ In Sec. 
Ill we apply the method of images to explain the vortex 
dynamics and we give some concluding remarks in Sec. 
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II. VORTEX CREATION AND DYNAMICS 
NUMERICALLY 

We assume initially a condensate trapped in a Gaus- 
sian optical trap. By shaking the trap rapidly we split the 
condensate into two spatially separated parts jl6| . First 
we slowly increase the amplitude of the periodic shak- 
ing along the chosen s-axis The effective potential 
evolves smoothly from the Gaussian well into a trap with 
two separated minima; this is the time-averaged poten- 
tial. The condensate follows adiabatically the changes 
of the time-averaged potential, splitting into two parts 
(Fig. |l|). By adding a similar process along the y-axis 
the time-averaged trapping potential evolves towards a 
torus. Thus the two separated condensate parts combine 
eventually into a torus. Before merging the two halves, 
we imprint on one of them (with a strongly detuned laser 
field 0) a change of phase of the wavefunction, A<f>. 
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At sufficiently low temperatures, the system is well de- 
scribed by the corresponding Gross-Pitaevskii equation 

ifi^- = -|-V 2 * + y* + NU Q \f\ 2 ty. (1) 
at 2m 

Here <J/ is the condensate wavefunction, V(x,y,t) is the 
trap potential JlTj (not the time-averaged one), N is the 
number of atoms and Uq = Anf^a/m, where a is the 
s-wave scattering length and m is the atomic mass. In 
all simulations we have assumed a sodium condensate 
(a = 2.75 nm) with the particle number N = 10 5 . The 
width of the trap in x- and y-directions was 10 fim and 
its size in the z-direction L z — 2 /jm. The trap depth was 
set to 4 fiK. In our simulations the duration of each shak- 
ing stage was typically 50 ms. The creation of the torus 
can be seen in the first steps of Fig. [p. With the chosen 
parameters the trap ground state energy in z-direction is 
considerably larger than the mean field energy nllo at the 
toroidal trap (n is the atomic density). This fact allows 
us to restrict ourselves to a two-dimensional geometry. 

The coherence time for the applied phase shift depends 
on the details of the shaking. If we express the wavefunc- 
tion as \E' = Rexp(i(f)) and assume the Thomas-Fermi 
limit, i.e., ignore terms due to kinetic energy in Eq. (Q), 
we obtain (assuming that R changes sufficiently slowly) 

d</> _ V(r,t)+NU R(r,t) 2 

Since in the Thomas-Fermi limit the right-hand side of 
Eq. (||) becomes independent of r, the phase evolution is 
identical everywhere. During the time between applying 
the phase shift and the merging into a torus, the phase 
has changed by a constant amount over the condensate 
halves and thus there will be a discontinuity in the phase 
across the merging points, corresponding to the applied 
phase. If, on the other hand, the condensate exhibits 
sloshing, the kinetic energy term makes the phase evo- 
lution r-dependent and the applied phase shift loses its 
meaning. Even a small asymmetry between the two parts 
is then enough to mask the applied phase shift on the 
time-scale of ~ 100 ms, which in effect leads to a total 
loss of the control of the vortex production. In our sim- 
ulations the density distribution followed the Thomas- 
Fermi limit closely and, therefore, the phase evolution 
was almost identical everywhere. 

When the phase-shifted condensate parts form a torus, 
the discontinuity in phase is quickly transformed into a 
dark soliton, since the system can adapt to the phase 
discontinuity only by reducing its density. In general, 
such a kink-like state in a trapped condensate is dy- 
namically stable if the mean-field energy nUo is roughly 
smaller than the characteristic trap energy in the direc- 
tion given by the soliton front Jl8| |, which in our case 
corresponds to the typical energy of the trap in the 
radial direction. With the parameters from our setup 
JlTj and n ~ 1.7 • 10 14 cm~ 3 , this condition is vio- 
lated. Under such conditions the soliton front under- 
goes snake instabilities which can eventually decay into 



vortex-antivortex pairs (or vortex rings in 3D configu- 
rations), see e.g. |10| , |l8| -[22| . This process has been ex- 
perimentally observed in non-linea r op tics [ p3| , as well 
as, very recently, in matter waves p3j . Since the front 
perturbation with a wavevector k = l/(i/2£), where 
£ = 1/ v / 47rna is the healing length, will grow fastest jl8) , 
the distance between the generated vortices is expected 
to be about d ~ ir/k. If k > l/£ the normal modes 
are dynamically stable, i.e., if the torus is narrower than 
about 7r£ the soliton is dynamically stable. 

The depth of the soliton depends on the size of the 
phase discontinuity. It must be large enough for vortex 
generation. Otherwise the resulting soliton can decay 
into other excitations, or even remain stable. For ex- 
ample, a 7r/2 phase shift did not lead to generation of 
vortices in our simulations. Due to the 2tt degeneracy of 
the phase, the phase shift tt corresponds to the largest 
possible discontinuity. But one should note that the di- 
rection of the phase gradient over the discontinuity affects 
the subsequent evolution. If the phase difference is equal 
to tt, the 2tt degeneracy makes the direction of the gra- 
dient ambiguous. In our simulations the motion of vor- 
tices was initially indeterminate and thus very sensitive 
to small disturbances, when A</> ~ tt. For A<fi = 0.9tt (as 
in Fig. Q), however, we obtain a much more predictable 
behavior. Due to symmetry, the mirror image of this 
situation is obtained with Acj) = I.Itt (the lower vor- 
tex moves clockwise in the beginning, and the upper one 
counterclockwise in Fig. |l|). 

In our simulations each soliton decays into a pair of 
vortices of opposite circulation. However, one of the vor- 
tices, located close to the edge of the condensate, has a 
very short lifetime, and decays within a few ms at the 
edge. The survival of only one vortex per soliton is not 
surprising, considering that d ~ I, where I is the width of 
the torus. For wider torii, i. e., longer solitons, our sim- 
ulations show a decay of each soliton into many vortex- 
antivortex pairs, as expected [ pp| . Similar behaviour has 
been reported in simulations of condensates in rectan- 
gular boxes (24|]. In our situation the condensate wave 
function remains symmetric with respect to reflection to 
the horizontal axis (see Fig. [l]). This symmetry leads to 
correlations between vortex production at the two phase 
discontinuities, located on opposite sides of the torus. 
We are always left with two vortices (one vortex per soli- 
ton) of opposite circulations. These vortices have also 
opposite vertical and equal horizontal components of the 
velocity, and they are therefore automatically set on a 
collision course. This reflects the conservation of angular 
momentum (we have verified that the shaking process or 
the adiabatic merging of the separated condensate parts 
do not bring any extra angular momentum into the sys- 
tem) 

Figure |l| shows the time evolution obtained from nu- 
merical simulations. The created vortices bounce from 
each other when they approach. This process repeats it- 
self if we continue simulations beyond t — 300 ms. In a 
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nondissipative homogeneous system two vortices of op- 
posite circulation separated by a distance greater than 
the healing length move parallel, since each vortex will 
move with the velocity of the other one 0. Clearly the 
vortex dynamics in a torus trap is more complex and 
completely different, ft can be understood with the help 
of the images method. 



III. IMAGES METHOD 

The images method has been recently successfully em- 
ployed to describe both a single vorte x |25| , and vortex 
arrays dynamics in circular 2D traps fl26[. We consider 
in the following a stiff torus (with infinitely high and 
stiff walls), i.e. we neglect the influence of the inhomoge- 
neous trapping potential on the vortex dynamics |p7| , p8f . 
Additionally, we assume that the distance between the 
vortices is greater than the size of the vortex cores. In 
absence of friction the vortex velocity equals the super- 
fluid velocity at the vortex position. Such a velocity 
field, which is induced by the presence of the vortices, 
has to fulfill the constraint of vanishing normal compo- 
nent at the torus boundaries. To fulfill this constraint 
it is enough to consider for each vortex an infinite num- 
ber of fictitious image vortices, each with appropriate 
position and circulation. Each vortex, either real or im- 
age, contributes to the supcrfluid velocity field vsf(t) by 
(«i/27r) z x (r — rj)/|r — | 2 , where and ft, refer to its 
position and circulation, respectively. 

Let us consider the case of a torus with an inner radius 
Ri , an outer radius R2 , and one vortex with circulation 



k located at r (in a frame centered at the origin of the 
torus). In order to cancel the normal component of the 
superfluid velocity on the boundaries, we have to con- 
sider two different families of image vortices. For the first 

family, an image vortex is placed at rj TO j = — - with cir- 
culation — k as well as one image vortex with circulation 
k at the origin pjJ. The image vortex at rj m i will induce 

a velocity component normal on the outer circle, so in ad- 

R 2 R 2 

dition we need an image vortex at Yi m 2 = 7-^-7 = ■j=j§r 
with circulation k. In turn this image will induce some 
normal component on the inner circle, so two more image 
vortices are required: one at Ti m ^ = 7-^^ = j^:^ with 
circulation — k the other one at the origin with circula- 
tion k and so on and so forth. For the second family, one 

considers an image vortex at r'; m i = with circula- 
tion —k. This image vortex gives a normal component 
of the velocity on the inner circle so we need in addition 

R 2 R 2 

one image vortex at r i m 2 — rr^- = -nir with circu- 
it 1 -"-2 

lation k as well as one image vortex at the origin with 
circulation — k. For this we require another image vortex 

at r'i m -i — -p 2 -- = -3?-- with circulation — k and so on. 

For the second real vortex in the torus we have to repeat 
the same formalism to get an image configuration, which 
ensures that the boundary condition is fulfilled. 

That way we can express the total superfluid velocity 
field in the torus (for arbitrary circulations), in particular 
the velocity at the locations of the real vortices as a func- 
tion of their mutual positions. Since the latter equals the 
vortex velocity we get the following first order equations 
for the vortex positions that can be easily computed 

\ 

+ Ri <-> R 2 

J 
\ 

+ Ri «-» R 2 

J 

e%» i,3 = l,2 (3) 




where Ri <-> R2 indicates a term equal to the previous 
one but interchanging i?i and i?2- The results obtained 
with the image method are presented in Fig. |^. As initial 
conditions we take two vortices located on opposite parts 
of the torus and displaced symmetrically from the cen- 
tral point between the torus boundaries. For the sake of 
clarity our initial conditions have been chosen such that 
the vortex trajectories are wider, i.e., they explore more 
space between the torus boundaries than in the numerical 



simulation from Fig. |lj. But with appropriate initial vor- 
tex positions one can achieve a nearly perfect agreement 
with simulations of Eq. (Q) concerning the shape of the 
trajectories, and the period of their oscillations (~ 300 
ms). 
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IV. CONCLUSIONS 



So far, we have assumed that the condensate is at zero 
temperature, i.e., we have used Gross-Pitaevskii equa- 
tion (|l|) to model the condensate, neglecting any dissipa- 
tion. The presence of non-condensate atoms is expected 
to result in dissipative Magnus forces and vortex decay. 
But experimental studies so far show that a realis- 
tic time scale for vortex decay is a few seconds, which 
is longer than the time scale of our scenario. We tested 
the stability of our system by adding a small, constant 
imaginary term into Eq. ([!]). It did not modify the re- 
ported behavior unless the magnitude of the term was too 
strong to be considered a perturbation. Nevertheless, the 
problem of the effects of dissipation in our model is very 
challenging. The presence of dissipation might cause the 
vortex-antivortex pair either to collide with the trapping 
boundaries, or to mutually annihilate. 

Another complication arises from the fact that cur- 
rently the spatial resolution of the measurement of the 
condensate density profile is limited. Usually this prob- 
lem is avoided by letting the condensate expand rapidly 
by removing the trapping potential. In a quasi-2D traps 
the optical depth will be an issue. Therefore some 
other detection scheme, such as matter wave interfer- 
ence J7|,p2|-p5[, may have to be considered. 

Summarizing, we have presented a method for creating 
vortices in a toroidal geometry in a controllable way. This 
method allows for the analysis of the interaction of vor- 
tices in such a geometry, which is significantly modified 
by the geometry of the trapping potential. The realm of 
possible phenomena that can be generated and controlled 
by merging the split condensate parts into a torus is not 
by any means limited by those discussed here, e.g. if we 
form the torus rapidly, other excitations are also created, 
which interact with the vortices and affect their dynam- 
ics. 

In the analytical estimation of the vortex orbits we 
have neglected for simplicity the effects of the inhomo- 
geneity of the trapping potential. Such inhomogcncity 
will result in an additional contribution to the velocity 
of the vortex in the direction perpendicular to the den- 
sity gradient, as discussed in Ref. p^] (a detailed analysis 
will be the subject of a later publication). On the other 
hand, our calculations have been constrained to quasi- 
2D traps. Very recently the issue of lower dimensional 
BEC has been subject of great interest. In particular 
2D (and even ID) condensates have been experimentally 
observed Therefore the situation considered in 

the present paper is experimentally justified and feasi- 
ble. However, the analysis of 3D geometries constitutes 
a challenging problem. We have made some short simu- 
lations to study the dynamics of vortex-antivortex (line) 
interactions in a 3D pipe geometry. The results we ob- 
tain are very similar to those presented here. In a proper 
3D geometry, the solitons decay into vortex rings, as ob- 
served very recently in the case of cylindrical traps p3| . 



The propagation, deformation and interaction of vortex 
rings will be the topic for future studies. 
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FIG. 1. Numerical simulation of the creation and dynamics of a vortex-antivortex pair in a Bose-Einstein condensate. The 
condensate is split into two parts, and then a phase shift A<f) — 0.9n is applied to the left part of the wavefunction. As the two 
parts merge to form a torus two solitons which subsequently decay into vortices are created. The two surviving vortices move 
along the torus and bounce from each other as they collide (minimum distance at t = 186 ms). 
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This figure "figl.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0005136v3 



